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SUMMARY 


A  free  flight  experiment  has  been  conducted  on  a  hemispherical 
shell  to  test  a  new  technique  for  deriving  outer  surface 
temperature  and  heat  transfer  distributions  from  inner  surface 
temperature  measurements  recorded  from  three  arrays  of 
thermocouples .  The  hemispherical  shell  which  was  mounted  on 
the  front  of  a  rocket  vehicle  reached  a  maximum  Mach  number  of 
5.5  at  4  km  altitude  after  7  s  of  flight.  Derived  heat  transfer 
profiles  are  characteristic  of  boundary  layer  flow  which  is 
laminar  in  the  stagnation  region,  changing  to  turbulent  within 
the  first  twenty  degrees  of  arc.  For  the  maximum  heating 
profile,  Stanton  numbers  varied  between  0.0037  at  30°  from  the 
stagnation  point  and  0.0025  at  60°.  The  maximum  free  stream 
Reynolds  number  based  on  body  diameter  was  26  x  10  .  Agreement  ^ 
with  laminar  stagnation  point  heat  transfer  theory  for  steady 
flow  is  good  except  for  just  after  maximum  heating  where  the 
experimental  results  are  some  twenty  five  per  cent  higher.  It 
is  concluded  that  the  new  technique  is  successful  but  that  the 
circumferential  temperature  and  heat  transfer  resolutions  could 
be  improved  by  increasing  the  number  of  thermocouples  in  each 
array  from  seven  to  at  least  nine. 
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1.  INTRODUCTION 

This  report  describes  the  results  of  a  free  flight  rocket  experiment  which, 
although  conducted  in  1971,  has  nevertheless  not  been  formally  recorded  until 
now.  This  experiment  marked  the  end  of  a  long  program  of  collaboration  between 
the  Royal  Aircraft  Establishment,  UK,  and  Weapons  Research  Establishment.  Under 
this  collaborative  program  a  hypersonic  research  rocket  was  developed  and  a  series 
of  these  rocket  vehicles  was  used  to  conduct  experiments  aimed  primarily  at 
investigating  heat  transfer  and  pressure  distribution  on  wings  and  noses  of  various 
shapes.  In  the  experiment  described  in  this  report  a  new  technique  has  been  used 
to  deduce  the  heat  transfer  to  the  outer  surface  of  a  hemispherical  shell-shaped 
nose  in  flight  from  measurements  of  its  inner  surface  temperature.  Of  particular 
interest  was  the  peak  turbulent  heating  and  the  stagnation  point  heating.  This 
shell  was  282  mm  in  diameter,  3.18  mm  thick,  and  made  from  Inconel  steel. 

To  deduce  heated  surface  conditions  from  inner  temperature  measurements  is  not 
altogether  straightforward  and  the  problem  is  a  member  of  the  class  known  as 
inverse  problems,  where  the  quantity  to  be  deduced  cannot  be  expressed  explicitly 
in  terms  of  the  measured  quantity.  The  greater  the  thickness  of  the  shell,  the 
greater  the  difficulty  becomes,  until  a  solution  eventually  becomes  inpossible. 
However,  during  other  work  carried  out  within  the  Weapons  Research  Establishment 
on  designing  and  analysing  heat  transfer  experiments  some  techniques  have  been 
discovered  which  delay  the  onset  of  these  difficulties  enabling  thicker  shells  to 
be  used.  These  techniques  have  been  put  to  use  in  the  present  experiment  where 
the  Inconel  steel  hemispherical  shell  was  thin  enough  to  make  an  inverse  problem 
solution  feasible  under  conditions  of  high  turbulent  heat  transfer,  but  thick 
enough  to  escape  melting  when  boosted  to  Mach  5.5  at  altitude  3.7  km  in  6.8  s. 

More  light  on  the  nature  of  the  difficulty  with  the  inverse  heating  problem 
as  such  can  be  shed  by  considering  one-dimensional  heat  flow  across  a  wall  heated 
on  one  surface  and  insulated  at  the  other.  Here,  using  the  Fourier  heat  conduction 
equation  as  the  mathematical  model,  leads  to  an  integral  equation  which  relates 
the  insulated  surface  temperature  explicitly  to  an  integral  containing  the  heated 
surface  temperature.  This  integral  equation  cannot  be  inverted  analytically  to 
give  the  heated  surface  temperature  explicitly  in  terms  of  the  insulated  surface 
temperature  which  would  be  the  measured  quantity  in  the  present  problem.  When  the 
heated  surface  temperature  variation  is  approximated  by  a  tabulated  function  with 
values  given  at  regular  time  intervals  of  8 ,  the  integral  equation  reduces  to 
algebraic  recurrence  relationships  which  can  be  inverted  to  give  heated  surface 
heat  flow  and  temperature  in  tabular  form  expressed  explicitly  in  terms  of  tabulated 
insulated  surface  measurements  and  previously  calculated  heated  surface  temperatures. 
One  might  now  think  that  this  approximate  solution  can  be  taken  to  the  accuracy 
required  by  making  the  time  step  8  small  enough.  This  turns  out  not  to  be  the 
case.  The  effect  of  making  6  smaller  is  to  introduce  higher  frequency  terms  in 
representing  the  inaccuracies  of  the  algebraic  process  and  insulated  surface 
temperatures  magnify  in  proportion  to  their  frequencies  when  translated  into 
heated  surface  temperatures.  This  magnification  is  dependent  not  only  on  the 
time  step  8  but  also  on  the  ability  of  the  wall  to  diffuse  heat,  the  diffusivity 
a,  and  the  wall  thickness  1.  The  relevant  nondimensional  frequency  parameter  for 
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tabulated  data  is  a  8/1  and  below  a  certain  critical  value  of  this  parameter  the 
algebraic  inversion  process  can  become  violently  unstable.  On  the  other  hand,  if 
the  time  step  8  has  to  be  made  too  large  to  avoid  this  instability,  the  data  will 
lose  its  sharpness  and  features  of  the  temperature  profiles  will  be  lost.  Fortunate¬ 
ly  there  is  a  way  around  this  problem  which  allows  the  time  step  8  to  be  chosen  sc 

that  the  temperature  profile  features  are  retained  and  the  onset  of  numerical 
instability  is  delayed.  This  is  brought  about  by  what  is  termed  the  wal 1 -splitting 

technique  which  uses  the  fact  that  for  one-dimensional  heat  flow  a  wall  can  be 

treated  exactly  as  an  equivalent  number  of  walls  of  lesser  thickness.  This  then 
has  the  effect  of  reducing  the  thickness  1  for  each  sub-wall,  thus  increasing  the 
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parameter  a  5/1  without  increasing  5.  Although  this  seems  to  have  the  disadvantage 
that  a  number  of  walls  have  to  be  treated  instead  of  just  one,  the  accumulative 
error  incurred  by  repeated  stable  algebraic  inversion  processes  does  remain  small. 
The  foregoing  discussion  relates  to  one -dimensional  heat  flow  and  the  next  question 
is  how  to  relate  this  to  the  assumed  two-dimensional  heat  flow  in  the  hemispherical 
shell. 

The  hemispherical  shell  is  assumed  to  be  heated  axi symmetrically  about  the 
longitudinal  axis  of  symmetry  normal  to  its  base.  Thus  the  mathematical  model  of 
the  heat  conduction  is  the  two-dimensional  Fourier  heat  conduction  equation  in 
spherical  polar  coordinates  where  heat  flows  radially  and  circumferentially  in 
meridian  planes  containing  the  longitudinal  axis  of  symmetry.  The  components  of 
temperature  arising  from  the  circumferential  heat  flow  are  separated  out  by  an 
expansion  of  the  temperature  solution  into  an  infinite  series  of  Legendre 
polynomials  whose  coefficients  are  functions  of  time  and  radial  position.  These 
coefficients  satisfy  an  equation  somewhat  similar  to  the  equation  satisfied  by 
temperature  in  one-dimensional  heat  flow  across  plane  walls.  Thus  a  somewhat 
modified  form  of  the  wall-splitting  technique  for  plane  walls  can  be  adapted  to 
these  coefficients  for  application  to  hemispherical  shells. 

So  far  the  question  of  noise  on  the  experimental  temperature  measurements  has 
not  been  mentioned.  Because  the  derivation  of  heated  surface  temperature  from 
insulated  surface  temperature  is  essentially  a  magnification  process  the  presence 
of  noise  on  the  insulated  surface  measurements  is  very  troublesome  and  raw  data 
must  be  smoothed  before  being  used  in  the  heated  surface  temperature  derivation 
process.  The  temperature  profiles  measured  in  this  experiment  are  very  well 
defined,  there  being  about  10  s  of  record  for  each  thermocouple  sampled  at 
80  points/s.  This  type  of  data  is  particularly  suited  to  overall  smoothing  by 
Fourier  analysis,  using  sine  series.  A  linear  temperature  variation  with  time  is 
combined  with  the  temperature  record  to  be  smoothed,  the  two  constants  being  chosen 
such  that  the  value  of  the  combined  function  is  zero  at  each  end.  When  this 
combined  function  is  regarded  as  an  odd  function  of  time,  both  it  and  its  first 
derivative  are  continuous  at  the  end  points  with  the  consequence  that  the  resulting 
Fourier  sine  series  expansion  tends  to  converge  as  n-3  where  n  is  the  term  number. 

At  least  this  is  true  for  the  smooth  part  of  the  function  but,  when  the  noise  begins 
to  dominate  in  the  higher  frequency  components,  this  tendency  to  converge  disappears 
and  the  higher  order  terms  assume  random  amplitudes  within  the  band  of  noise  level. 
This  distinctive  behaviour  enables  noise  to  be  separated  from  the  smooth  part  of 
the  function  in  Fourier  smoothing.  In  this  experiment,  with  about  800  points  per 
record  there  is  enough  redundant  information  to  be  able  to  least-squares  fit  a  300 
term  Fourier  sine  series  to  each  record.  Another  advantage  of  Fourier  smoothing 
in  free  flight  experiments  is  that  it  enables  spurious  frequencies  to  be  detected 
such  as  those  which  may  arise  from  rotary  switches  or  structural  vibrations. 

The  object  of  this  experiment  is  to  provide  a  practical  test  for  the  techniques 
previously  mentioned  and  at  the  same  time  to  acquire  some  meaningful  information 
on  the  surface  heat  flux  and  temperature  distributions  around  a  hemispherical  nose 
during  low  altitude  supersonic  flight.  In  the  next  section  a  description  of  the 
rocket  test  vehicle,  its  performance  and  instrumentation,  is  given.  A  fuller 
account  of  the  theory  and  method  together  with  references  is  given  in  Section  3 
and  a  discussion  of  the  results  and  an  assessment  of  the  techniques  used  is  given 
in  Section  4. 


2.  ROCKET  TEST  VEHICLE  AND  INSTRUMENTATION 

The  hemispherical  shell  tested  formed  the  nose  of  the  rocket  vehicle  shown  in 
figures  1  and  2.  Further  details  of  this  rocket  vehicle  together  with  reasons  for 
using  it  are  given  in  reference  1.  Trajectory  data  were  obtained  from  radio  Doppler, 
Cont raves  kinetheodolites ,  and  radar.  Vinten  cameras  were  also  used  to  monitor  the 
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behaviour  of  the  vehicle  itself  during  flight.  The  vehicle  performed  satisfactorily 
and  at  apogee  the  hemisphere  attached  to  the  instrumentation  compartment  was 
separated  from  the  boost  motor  to  be  recovered  intact  on  the  ground  some  21  km 
downrange  after  deploying  first  a  drogue  and  then  a  parachute.  The  sequence  of 
events  leading  to  recovery  is  shown  in  figure  3.  Figure  4  shows  Mach  number,  free 
stream  Reynolds  number,  altitude,  and  stagnation  temperature  histories  for  the 
first  nine  seconds  of  flight.  Figure  S  shows  a  photograph,  taken  shortly  after 
impact  on  the  Woomera  range,  of  the  hemisphere  and  instrumentation  compartment 
with  the  recovery  parachute  still  attached.  No  discernible  damage  was  suffered 
by  the  hemisphere  which  with  thermocouples  intact  is  avilable  for  further 
experiment,  either  wind  tunnel  or  flight. 

The  diameter  of  the  hemispherical  shell  was  282  mm,  its  thickness  was  3.18  mm 
and  it  was  made  out  of  'Inconel  600,'  a  material  chosen  because  its  thermal 
diffusivity  varied  less  with  temperature  than  other  available  materials  and  the 
mathematical  model  used  for  the  heating  analysis  was  based  on  constant  diffusivity. 

To  ensure  the  presence  of  fully  turbulent  flow  on  at  least  part  of  the  hemisphere 
the  outer  surface  was  finished  to  a  roughness,  nowhere  less  than  0,64  Pm.  Nineteen 
chromel-alumel  thermocouples  were  spot  welded  to  the  inner  surface  of  the  shell 
along  three  generators  in  planes  120  degrees  apart,  their  positions  being  shown  in 
figure  6  together  with  four  pressure  monitoring  holes  connected  to  a  differential 
pressure  yawmeter.  To  reduce  heat  conduction  losses,  a  durestos  insulating  ring 
one  inch  thick  was  inserted  between  the  base  of  the  hemisphere  and  the 
instrumentation  compartment  and  a  hemispherical  radiation  reflecting  surface  was 
fitted  within  the  shell.  Again,  the  mathematical  model  of  heat  conduction  in  the 
shell  assumed  that  no  heat  crossed  the  inner  surface  or  the  base. 

The  instrumentation  compartment  carried  a  24  channel  465  MHz  telemetry  system 
which  transmitted  data  to  ground  receivers.  Nineteen  channels  were  used  to  transmit 
temperature  data  from  each  thermocouple  at  the  rate  of  eighty  points  per  second, 
one  channel  was  used  for  the  synchronising  pulses  needed  to  identify  data  on  records, 
another  for  inductance  calibrations,  and  the  remaining  three  were  sub -commutated  so 
that  each  could  sample  twenty  four  different  inputs.  One  of  the  sub -commutated 
channels  carried  data  from  all  the  inductance  transducers  used  to  measure  lateral 
and  longitudinal  acceleration  and  yaw  meter  differential  pressures.  The  other  two 
sub-commutated  channels  carried  data  from  two  sun  cells,  from  the  recovery  system 
monitor,  and  voltage  calibrations  including  inputs  from  a  cold  junction  and  ten 
thermistors  for  subsequent  calibration  of  thermocouple  output.  Telemetered  data 
after  being  received  on  the  ground  was  recorded  in  analogue  form  on  magnetic  tape 
and  subsequently  converted  to  digital  form  for  input  to  a  computer  program  which 
used  the  in-flight  calibrations  to  produce  time  histories  of  the  physical  quantities 
being  measured. 

No  rolling  motion  was  evident  during  the  first  nine-seconds  of  flight  from 
either  the  sun  cell  records  or  from  the  Vinten  camera  film.  Incidence  history  as 
deduced  from  the  differential  pressure  yawmeter  was  obviously  wrong;  the  results 
showed  an  incidence  of  half  a  degree  for  the  first  five  seconds  followed  by  a 
steady  increase  of  about  one  degree  per  second  which  would  have  led  to  break  up  of 
the  vehicle's  structure.  The  raw  temperature  measurements  from  each  of  the  nineteen 
thermocouples  are  shown  in  figures  7  to  14  covering  the  first  nine  seconds  of  flight. 


3.  THEORY  AND  METHOD 

The  mathematical  model  of  the  heat  conduction  used  in  this  analysis  is  that  of 
a  hemispherical  shell  of  inner  radius  a  and  outer  radius  b,  having  constant  thermal 
diffusivity  a.  The  symmetry  of  the  external  heating  is  assumed  to  be  such  that 
heat  can  flow  only  in  the  r  and  6  directions  shown  in  figure  15;  then  the  temperat¬ 
ure  at  a  point  r,  at  time  t,  T  (M,  r,  t),  where  M  =  cos  6 ,  is  given  by  the  equation 
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f^)*  h  It 

Details  of  the  derivation  of  this  equation  can  be  found  in  reference  2.  It  is 
well  known  that  solutions  of  this  equation  can  be  expanded  in  series  of  Lengendre 
polynomials  in  M  where  the  coefficients  are  functions  of  r  and  t.  A  convenient 
way  of  arriving  at  this  result  for  the  purpose  of  the  present  analysis  is  to 
introduce  the  integral  transform 


1 

^2n  (r,t)  =  ^  T  CAi,  x,  t)  P2n  (M)  d  »  (2) 

0 

In  equation  (2)  the  P2n00  are  LeSen<*re  polynomials  of  order  2n  in  even  powers 

of  M  and  the  T2n(r,t)  can  be  considered  as  Legendre  components  of  temperature 

T  (M,r,t).  In  fact  temperatures  are  reconstituted  from  their  Legendre  components 
by  the  inverse  transformation 


T  (M ,  r ,  t ) 


OO 


(4n+ 1) 


2n 


00  T2n(r,t) 


(3) 


To  find  the  T2^(r,t),  the  integral  transform  given  by  equation  (2),  when  applied 
to  equation  (1)  gives 


1_ 

2 

r 


2n(2n+l)  T 

2  ^ 
r 


1_ 

a 


at 


(4) 


where  it  has  been  assumed  that  no  heat  passes  through  the  plane  8  =  which 

a  t  ^ 

implies  the  condition  =  0  at  M  =  0.  To  provide  for  this  condition  in  the 

experiment  a  durestos  ring  was  used  to  insulate  the  base  of  the  hemisphere  from 

the  adjoining  instrumentation  compartment.  It  should  be  noted  here  too,  that  T 

in  equations  (2)  and  (3)  can  be  replaced  by  the  radial  heat  flow  Qr  which  is  given 

9T  ~~  0T 

by  -K  —  (M,r,t)  with  Legendre  component  Q2n Cr> defined  by  -K  ^  (r,t),  where 

K  is  the  coefficient  of  thermal  conductivity  for  the  shell  material. 

In  this  experiment,  temperature  histories  were  measured  by  three  rows  of 
thermocouples  attached  to  the  inner  surface  r=a  of  the  hemisphere  and  are  shown 
in  figures  7  to  14.  Hence  in  order  to  use  a  one -dimensional  analysis  based  on 
equation  (4)  to  deduce  temperature  and  radial  heat  flow  at  the  outer  surface  r=b, 

it  is  necessary  to  first  find  the  Legendre  components  T2  (a,t)  of  the  measured 

temperatures.  This  was  accomplished  by  fitting  piecewise  quadratics  in  p  to  the 
time-smoothed  measurements  T  (M,a,t)  through  each  three  successive  M  positions 
along  each  row  of  thermocouples.  These  quadratic  expressions  for  T  (P , a , t )  were 

then  inserted  into  equation  (2)  from  which  values  of  UP  t0  n=6  were 

calculated.  Thus  only  a  truncated  version  of  equation  (3)  was  used  in  reconstitut¬ 
ing  temperatures  and  heat  flows  from  their  Legendre  components. 

The  method  used  for  deriving  ^  (b,t)  ^2n^,t:)  ^rom  ecluation  (4)  follows 


5 
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similar  lines  to  that  used  for  deriving  outer  surface  temperature  and  heat  flow 
from  inner  surface  temperature  history  measurements  taken  on  flat  plates.  The 
solution  of  equation  (4)  when  no  heat  passes  across  the  surface  r=a  can  be 
extracted  from  reference  (2)  or  obtained  directly  from  reference  (4),  and  is 


T2„l 


r,t) 


+  / b\ 
\  r/ 


(a»2n+)s  (r\2n+!s  (  1_\ 

\tI  +  \  a/  \1  +  2n/ 

/a\2n+Ji  *(b\2n+li  {  n 

\b/  \a/  \1  +  2nj 


T2n(b.t) 


S2n+h  (gs,r) 


3T2n(b,n  a(t-?) 


i 9? 

s  u 


d  f 


where 


{' 


B2nA  ^S’r)  =  J2n+/S  r)  \H  Y2n+/S  a)  ’  °s  a  Y'2n^  (°s  a) 


) 


(5) 


Y2n+j/s  r)  1  4  J2n+i/S  a)  “  °S  a  J' 2n*h  (°S  a) 


(6) 


J-  i  and  ,  are  Spherical  Bessel  Functions  of  the  first  and  second  kinds 
2n+^  2n+%  r 

respectively  and  s  is  the  sth  zero  of  B  1  (ffjb).  The  properties  of  Spherical 

Bessel  Functions  are  given  in  reference  3  where  they  are  shown  to  be  expressible 
in  the  form 

•>*,♦*«>  '  ^{f2n  (8)  Si"  1  '  f-2n-l(Z>  Cos  Z  } 

*  -^1|{f2n  ‘Z>  C0S  Z  *  f-2„-l<Z>  Si"  Z 

and  the  fm(Z)  are  finite  polynomials  in  negative  powers  of  Z  such  that 
fm-l(Z)  +  fm+1 (Z)  =  (2m+l)  Z_1  fm(Z) ,  (m=0,  ±1,  ±2,  ...  ) 
fO(Z)  =  Z'1  and  fl(Z)  =  Z-2 


Thus  the  Spherical  Bessel  Functions  for  the  hemispherical  shell  are  simply  related 
to  the  circular  functions  cos  Z,  sin  Z,  which  occur  without  the  fm  polynomials  in 
one-dimensional  flat  plate  heat  conduction  theory.  In  fact,  equation  (5)  is 
similar  in  form  to  the  flat  plate  solution  which  links  the  inner  insulated 
surface  temperature  with  the  outer  heated  surface  temperature. 


3.1  Derivation  of  Outer  Surface  Temperature 

The  difficulties  in  deriving  T2n(b,t)  from  through  equation  (5) 

with  r=a,  are  similar  in  nature  to  the  flat  plate  inverse  problem,  so  that  the 
devices  used  to  conquer  that  problem  can  also  be  applied  to  the  shell. 
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Because  there  is  no  explicit  way  in  which  T(b,t)  can  be  expressed  as  a 

continuous  analytic  function  of  T(a,t),  a  finite  difference  approach  is 

adopted  in  determining  T(b,t)  from  T(a,t)  through  equation  (5)  with  r=a. 

This  method  is  explained  in  detail  in  reference  (4).  In  short,  T^Ct^t)  is 

considered  to  consist  of  a  series  of  chords  joining  the  points  ^2n^,r 

to  T2n(b,r  5),  where  r  ranges  in  integer  values  from  1  to  m  and  m  5  =  t. 

Substituting  these  values  for  into  equation  (5)  gives  the  following 

result  for  the  Legendre  components  of  the  shell  inner  insulated  surface 
temperature 

/ b  \  2n+1 

T2n(a,t)  .(bf  (4-1)11) 

\a/  <«♦» 


TTn  s  3a  2n+V  s 
s-u  s 


m  r 


B2n+%  ^s,a) _ \  l_  fi  T 

JLL 


J.r'-lS^ 


o~a{  nfi-T) 


=1  r-1  6 


which,  on  evaluating  the  integral,  can  be  written  in  the  form 


m-1 

T2n(a,.5)  -jj 
__  r=l 


T  (b,m  6)  =  1  T  (a,m5)  -\  T.  (b,r5)  0 
2nk  1  J  7$  2n^  *  /  2nk  *  '  m-r+l,n 


where 


©i  =  (4n 
l,n  — 


/  b  \ 
+  1)  \  a  / 


■  1  7raS 


2n+  (2n+l)[ 


oo 

(I)  G 


o2a5 
(1-  e  s  ; 


a3  9  B 

^  s  2n+^(os,b) 


m-r+1 ,n 


(3'i> 


2  2  0 

a  aS  a ^  a  6 (m-r- 1) 
-  e _ )  e _ 

s3  J_  (a  ,b) 


and  the  are  the  zeros  of  equation  (6)  with  r=b. 


Thus  equations  (7),  (8)  and  (9),  allow  T2n  (b,m  5}  to  be  calculated  from  a 
series  of  algebraic  recurrence  relationships  involving  T2n  (a,m  5)  and  the 
previously  calculated  values  of  T^n  (b,r  5). 
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However,  the  calculation  of  an  outer  heated  surface  temperature  from  a 
known  inner  insulated  surface  is  essentially  a  magnification  process.  Thus 
any  noise  on  experimental  data  or  errors  inherent  in  numerical  procedures 
tend  to  become  magnified.  It  is  shown  in  reference  (2)  that  the  amplitude 
magnification  is  greater  on  the  higher  frequency  components  of  the  inner 
insulated  surface  temperature.  In  this  experiment  the  higher  frequency 
components  stem  from  noise  in  the  measured  temperatures  and  from  the  step 
size  8  used  in  the  numerical  calculation  procedure.  To  a  large  extent  the 
high  frequency  noise  in  the  temperature  measurements  has  been  removed  by 
a  Fourier  smoothing  process  which  fits  a  Fourier  sine  series  about  a  straight 
line  to  represent  each  temperature-time  record  and  truncates  the  series  at 
some  pre-assigned  cut-off  frequency.  Under  these  conditions  the  Fourier 

-3 

representation  of  a  continuous  signal  tends  to  converge  as  n  whereas  a 
discontinuous  noise  component  has  no  tendency  to  converge.  This  effect  is 
well  illustrated  in  figures  16,  17  and  18,  which  show  the  relationship  between 
amplitude  and  frequency  for  the  temperature  measurements  taken  at  0  =  120°, 

6  =  20°,  30°,  and  40°.  In  these  figures,  the  amplitude  noise  band  has  been 
chosen  to  include  amplitude  values  which  obviously  exhibit  random  behaviour. 
Starting  with  the  lower  frequencies,  the  standard  deviation  of  these  random 
amplitudes  approaches  a  limit  as  the  number  of  values  included  is  increased. 

A  cut-off  amplitude  corresponding  to  a  fixed  multiple  of  this  limiting  value 
of  the  standard  deviation  of  the  random  amplitudes  is  then  selected  and  only 
amplitudes  above  this  cut-off  value  are  retained  to  represent  smoothed 
temperatures.  For  this  experiment,  only  amplitudes  greater  than  six 
standard  deviations  of  the  random  amplitudes  have  been  retained.  The  smoothed 
temperature  was  represented  by  a  straight  line  plus  the  first  four  Fourier 
terms  for  0  =  20°,  whereas  the  first  six  terms  were  used  in  representing 
0  =  30°  and  0  =  40°.  These  smoothed  temperature  representations  are  shown 
superimposed  on  the  raw  data  in  figure  19.  Smoothed  temperatures  for  0=0, 

0  =  50°,  70°  and  90°  are  shown  in  figure  20  where  six  Fourier  terms  were  used 
for  0  =  50°  and  four  for  each  of  0  =  70°  and  0  =  90°.  Having  removed  the  high 
frequency  noise  from  the  raw  data  the  Legendre  components  of  the  smoothed 
temperatures  are  represented  as  tabulated  functions  with  values  given  at 
regular  time  intervals  8.  These  values,  when  introduced  into  equation  (7), 
which  is  the  algebraic  inversion  process  for  calculating  the  outer  surface 

temperature  components  T^Cb.m  have  the  effect  of  introducing  an  error  of 

frequency  8.  The  larger  the  thermal  thickness  of  the  shell  and  the  smaller 
the  time  step  8  the  greater  the  amplification  of  this  error  in  the  derived 
temperature  components.  This  behaviour  is  just  contrary  to  what  is  required; 
the  shell  needs  to  be  thermally  thick  enough  to  survive  the  high  turbulent 
low-altitude  heating  while  the  larger  the  time  step  8  the  coarser  the  numerical 
temperature  representation.  This  problem  can  be  eased  by  splitting  the  shell 
thickness  into  a  number  of  smaller  thicknesses  for  the  purpose  of  calculation. 
The  way  in  which  this  works  is  best  illustrated  by  briefly  describing  the  wall 
splitting  technique  as  given  in  reference  5  for  flat  plates.  Applying  the 


operational  transform  T  (x,p) 


heat  conduction  equation 


d  T(x,t) 
a  2 


(x,t) 


il(x,  t) 
dt 


dt  to  the  one-dimensional 


gives  the  solution 


T(l.p) 


_  T(o,p) 
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when  heat  is  received  at  the  surface  x=0  and  no  heat  passes  across  x=l . 
Equation  (10)  can  be  re-written  as 


T  (o,p)  =  T  (l,p)  Cosh  y/%  1 


where  Cosh  a  1  can  be  regarded  as  the  operator  which  when  applied  to  the 
transform  of  the  insulated  surface  temperature  gives  the  transform  of  the 
heated  surface  temperature*  Now  Cosh  >/j^  1  can  be  expanded  in  terms  of 
multiple  arguments  to  generate  operators  representing  walls  of  thickness  less 
than  1.  The  simplest  expansion  of  this  type  is  to  put  Cosh  \/^  1=2  Cosh^ 

1,  which  when  substituted  into  equation  (10)  gives 


v  a  2 


T  (o,p)  =  2  T  (l,p)  Cosh2/|l.  .  t  (1’P) 


This  result  is  interpreted  to  mean  that  the  temperature  T  (l,t)  is  applied 
to  the  insulated  surface  of  a  wall  of  half  the  original  thickness  and  the 
outer  surface  temperature  so  obtained  is  then  itself  applied  to  the  insulated 
surface  of  the  same  half-thickness  wall.  The  outer  surface  temperature 
obtained  from  this  second  operation  is  then  doubled  and  the  original  T  (l,t) 
subtracted  to  give  the  outer  surface  temperature  sought  after.  Because,  as 
previously  mentioned,  the  Spherical  Bessel  Functions  of  the  solutions  to  heat 
conduction  in  a  hemisphere  are  simply  related  to  the  circular  functions  of  the 
flat  plate  solutions,  the  wall-splitting  technique  can  be  applied  to  a 
hemispherical  shell  and  this  is  shown  in  reference  6.  In  fact,  the  computer 
program  described  in  reference  6  for  calculating  shell  outer  surface  temperatures 
from  inner  surface  temperatures,  as  used  in  this  experiment,  allows  the  shell 
thickness  to  be  subdivided  up  to  six  times. 

When  the  error  amplification  effects  become  large,  the  numerical  calculation 
of  outer  surface  temperature  by  a  process  such  as  equation  (7)  can  suddenly 
become  violently  unstable  causing  the  calculated  results  to  oscillate  wildly. 

This  feature  renders  the  onset  of  instability  very  easy  to  detect  in  running 
a  computer  program,  so  that  a  method  of  automatically  introducing  a  wall-splitting 
is  readily  arranged.  The  tendency  to  amplify  errors  manifests  itself  through 

©.  in  equation  (7)  where  it  is  found  that  ©„  decreases  as  — —  n  decreases. 

1-n  (b-a)2 

It  is  of  interest  to  note  that  both  Pope  and  Heron,  in  references  (7)  and  (8), 
found  independently,  that  to  avoid  instability  in  this  method  of  calculating 
outer  heated  surface  temperature  from  inner  insulated  surface  temperature  for 

flat  plates  of  thickness  1,  must  be  greater  than  2/3. 

r 

3.2  Derivation  of  Outer  Surface  Radial  Heat  Transfer 

The  Legendre  components  of  the  outer  surface  heat  transfer  are  determined 
from  the  outer  surface  temperature  components  as  calculated  by  equation  (7). 

This  is  achieved  by  putting 

Q2n  (b.t) 

which,  from  equation  (5),  gives 


9T 

9r 


(r,t) 


r=b 


r 


'***^W* 
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Q2n(b,t)  =  2n  (2n+l) 

IT  b 


3r  B2n+Js 


(b\4n+1  1  _ 

2n+l)  (111  -  1  >  T?n  Cb.t) 

\  /h\4n+l  ( 

£  (2n+l)  ( ^-j  +2n  j 

,r)]  r=b  fa 2n  Cb,n  e'asa  ^ 


(»  .b) 


^7q  **  as  To  B2n+^  **as 


Although  it  is  not  apparent  from  this  form  of  the  solution,  an  alternative 
solution  suitable  for  small  values  of  time  shows  that  the 

coefficient  of  _ 2n  in  equation  (11)  behaves  as  C  (t  -  ^  as  ?  ♦  t,  where 

9?  ax 

C  is  a  constant.  Thus  to  employ  the  approximation  that  2n(b,f)  can  be 

1  r  ^  *  _  3f 

taken  as  5"  )^2n(b,r  6)  "  ^2n(b,ai6)  f  for  r-l«  <f<r8 

in  equation  (11),  it  is  advisable  to  integrate  out  the  singularity  that  occurs 
as  ?-►  t.  This  is  achieved  by  integrating  both  sides  of  equation  (11)  from 
t  -  6  to  t,  giving 


^2n^b’^  df  =  2n  (2n-» 
K  b 


Jr . , 


*h  (°s’r)] 


(2n+l) 


t  r 


T2n(b,f)  d  f 


dr  2n 


h  o  _  .  (o  ,b) 

— f  s  3a  2n+b  s*  J 
s=0  s 


3T,  (b,n  e^sa$V 


d  ?  jd  f 


t-5  0 


After  changing  the  order  of  integration  in  the  double  integral,  one  integration 
can  be  performed  and  when  dT^  ^  j.  ^  is  replaced  by 

5T7  1 


T  (b,rS)  -  T  (b,r-l  6  >  for  r-1  8<  f<  r  6  and  /  Q  (b,f )  d  f  and 


T2n(b>f)  d  f  are  approximated  by  j  s  Q2n  Cb > + 
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A 


T2n(b,t)  +  T,n  (b,t-  6)  (  respectively,  equation  (12)  becomes 


2n 


Q2n(b-"{)  ’  <!2n(b,m"1  S>  *  T2n(b,»S)  Ml  n  *  T2n(b,m-1  «)  M2  n 
m-2 


•D 


L  (b,r8)  M  . 

2n  m-r+1,  n. 


r=  1 


where 


vn  \v  p 


3r  ^2n+%  s 


.5  3 


■r>Jr.b  [  e'°s  “*-l  .  .8  J 


s=0  °s  da  S2n+h  *'<?s’'3^ 


ter1. 

>(! 


+  2n  (2n+l)  ~ 


4n+l 


M 


4K  V't*  B2  , 
=  272  )  ldr  2n+l 

“  5  L.S  L.  » 


2n+^  (Qs 


•*»] ~ 


'o~  r—  B.  ,  (O  ,b) 
s=0  s  do  2n+h  s’ 
s 


+  2n  <2n+1>  F 


(2n+l)|  d  +2n 


2  2 

k  i  ^  -  -a  a  6  -2a  a  5  2  * 

b  /  2-3e  s  +  e  s  -  a  a  6 

s 


ter1., 


,.\4n+l 

(2n+D(|-)  +2n 


and 


V-...  ■  ^ffewvj  -  i':1 

^  L_  B-wu  (<>.*>) 


(13) 


(14) 


(15) 


(16) 


for  m-r+  1  >  2 . 


The  evaluation  of  ^2n+h^° s*T^  *  its  first  derivatives  with  respect  to  r  and 

a^,  and  its  zeros,  is  given  in  detail  in  references  (4)  and  (9).  In  reference 

(6)  a  Fortran  II  computer  program  is  described  for  the  inverse  problem  of  the 
heated  shell  using  the  wall-splitting  technique.  Reference  (9)  describes 
another  computer  program  for  calculating  the  surface  temperature  and  heat 
input  to  a  hemispherical  shell  from  a  prescribed  outer  surface  variation  of 
heat  transfer  coefficient  and  recovery  temperature. 


i 
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4.  RESULTS  AND  DISCUSSION 

The  raw  temperature  data  measured  by  each  of  the  nineteen  thermocouples  is 
shown  in  figures  7  to  14.  All  the  data  except  that  measured  at  the  three 
stations  for  0=20°  is  compatible  with  axisymmetric  flow  around  the  hemisphere. 
When,  as  is  shown  in  figure  11,  the  temperature  measurements  for  0=0,  0=20°,  are 
superimposed  on  those  for  0=120°,  0=20°,  a  marked  difference  is  apparent.  The 
temperatures  at  0=0  rise  more  steeply  than  those  for  0=120°.  The  incidence 
required  to  account  for  such  a  temperature  difference  would  be  more  than  the 
rocket  vehicle  could  withstand  and  the  thermocouple  at  0=0,  0=20°,  was  checked  in 
the  laboratory  after  the  trial  and  found  to  be  in  perfect  working  order.  This 
then  leaves  boundary  layer  transition  as  an  explanation.  Earlier  boundary  layer 
transition  in  the  0=0  meridian  plane  than  in  the  0=  ±120°  planes  due  perhaps  to 
different  surface  roughness  distributions  is  in  keeping  with  observed  boundary 
layer  behaviour  on  hemispherical  noses.  There  is  also  some  discernible  difference 
between  the  temperature  measurements  for  0=120°  and  0=  -120°  at  0=20°  and  although 
this  difference  is  not  nearly  so  marked  it  does  lend  weight  to  the  assumption 
that  asymmetric  boundary  layer  transition  occurred  in  the  neighbourhood  of  0=20°. 

For  the  purpose  of  smoothing,  each  record  of  measured  temperatures  was 
represented  by  a  Fourier  sine  series  fitted  about  a  stright  line.  Plots  of  the 
logarithm  of  amplitude  against  frequency  are  shown  for  the  Fourier  terms 
representing  the  temperatures  measured  at  0=120°,  0=20°,  30°,  40°  in  figures  16, 

17  and  18.  Similar  plots  were  obtained  for  all  the  other  measurements  and  they 
serve  to  show  the  range  of  amplitude  values  dominated  by  noise  and  the  range  which 
is  predominantly  signal.  From  the  amplitudes  representing  the  noise,  standard 
deviations  for  the  noise  in  the  Fourier  coefficients  have  been  calculated  and 
cut-off  frequencies  for  smoothing  were  chosen  such  that  only  those  frequencies 
having  amplitudes  which  exceeded  six  standard  deviations  were  retained.  As  a 
result  the  standard  deviation  of  the  error  due  to  noise  in  the  raw  data  among 
the  nineteen  records  varied  between  5.5°K  and  8.6°K,  whereas  the  standard 
deviation  of  the  error  in  the  smoothed  data  varied  between  0.49°K  and  0.82°K  when 
Fourier  terms  with  amplitudes  less  than  0.52°K  were  rejected.  Bias  errors  were 
found  to  be  quite  small  having  a  maximum  value  of  order  0.5°K.  Two  examples  of 
curves  of  Fourier  smoothed  temperatures  superimposed  on  raw  data  are  shown  in 
figures  19  and  20.  Most  of  the  noise  seems  to  occur  in  the  last  two  seconds  of 
each  record  when  ablating  silicon  rubber  used  as  thermal  protection  could  have 
passed  on  to  the  surface  of  the  telemetry  transmitting  aerials. 

<n - - _ i  _ _ _ _  _  _ i. _ i  i _  n-6 


The  smoothed  temperatures  were  represented  by 


2  V*'')  P2„(CoS«). 

n=0 


the  first  seven  terms  of  their  Legendre  polynomial  expansions.  Just  how  well  these 
first  seven  terms  represent  the  smoothed  temperatures  can  be  seen  from  figures  19 
and  20.  The  worst  fit  occurred  for  0=120°,  0=30°,  as  shown  in  figure  19  and  the 
result  for  0=  -120°,  0=30°,  was  much  the  same.  However,  the  representation  for 
all  the  other  data  displayed  the  same  quality  of  fit  as  shown  for  0=0,  0=50°,  70°, 
90°,  in  figure  20,  including  0=0,  0=30°.  This  raises  the  question  of  just  how 
many  terms  are  needed  in  the  Legendre  polynomial  expansions.  Certainly  the  seven 
terms  used  seem  adequate  for  all  but  the  0=  ±120°,  0=30°  cases  where  it  might  be 
expected  that  adding  about  two  more  terms  would  greatly  improve  the  representation. 
Where  a  fairly  large  number  of  terms  is  needed,  this  suggests  that  instead  of 
truncating  the  infinite  series  it  would  be  worthwhile  to  assume  a  finite  expansion 
in  the  first  place  and  find  the  Legendre  coefficients  from  a  system  of  simultaneous 
algebraic  equations  such  that  the  smoothed  thermocouple  readings  are  regenerated 
exactly.  This  process  would  be  much  simpler  than  calculating  the  exact  Legendre 
coefficients  through  their  integral  representations  as  defined  for  example  by 
equation  (2).  The  disadvantage  would  be  that  the  Legendre  coefficients  found  from 
the  finite  expansions  would  not  be  exact  but  they  do  have  the  useful  property  that 
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their  error  diminishes  as  the  number  of  terms  in  the  finite  expansion  increases. 

Even  in  attempting  to  find  the  exact  coefficients  there  will  be  errors  arising 
from  the  integration  of  products  of  cos11  0  and  a  function  represented  numerically 
or  approximated  by  a  series  of  suitably  simple  analytic  functions.  Before 
performing  another  heating  experiment  in  this  way  on  a  hemisphere  it  would  be 
worthwhile  carrying  out  an  accuracy  check  on  the  orders  of  finite  Legendre 
polynomial  expansions  required  to  represent  analytical  descriptions  of  typical 
temperature  profiles  tojwithin  given  accuracies.  ^ 

From  the  values  of  T2n(a,t)  the  Legendre  coefficients,  T2n(b,t),  of  the  outer 

surface  temperature  were  calculated  through  equation  (7),  enploying  the  wall¬ 
splitting  technique  where  necessary,  and  the  reconstituted  outer  surface 
temperatures  given  by 

n=6 

T  (cos  0,  b,t)  =  2  T  (b,t)  P2n(cos  0 ) 

n=0  "  A 

are  shown  plotted  in  figures  21,  22  and  23.  Each  row  of  thermocouples  along 
0=  ±120°  and  0=0  was  treated  independently  and  no  attempt  has  been  made  to  use 
the  fact  that  the  stagnation  point  thermocouple  was  common  to  each  array  and 
consequently  the  differences  in  the  derived  values  of  T  (b,t)  at  0=0  for  the 
three  values  of  0  serve  as  an  indication  of  the  error  in  the  method.  The  agreement 
between  T  (b,t)  at  the  stagnation  point  for  0=120°  and  0=  -120°  is  quite  good 
whereas  the  agreement  of  these  with  0=0  is  not  so  good.  As  shown  in  figure  11  the 
raw  temperature  data  profile  for  0=0  rises  more  steeply  than  that  for  0=120°  at 
0=20°  and  the  different  profile  has  produced  different  values  for  the  Legendre 
coefficients  of  which  only  seven  have  been  used  to  reconstitute  temperature.  This 
is  where  the  alternative  finite  Legendre  expansion  proposed  in  the  previous 
paragraph  would  have  ensured  a  value  of  T  (a,t)  at  the  stagnation  point  common  to 
each  row.  The  steeper  rise  in  the  raw  temperature  data  for  0=0,  0=20°,  is 
reflected  in  the  outer  surface  temperature  profiles  of  figure  21  when  compared 
against  the  profiles  for  0=120°  in  figures  22  and  23.  Another  feature  is  the 
oscillation  that  occurs  in  some  of  these  temperature  profiles  after  the  temperatures 
pass  through  their  maximum  values.  This  oscillation  is  undoubtedly  a  result  of  the 
Legendre  representation  by  seven  terms  where  the  last  term  retained  includes  a 
cos  0  to  the  twelfth  power  which  therefore  has  a  cos  12  0  component  with  a  period 
of  30°.  The  oscillation  does  have  a  period  of  roughly  30°  and  manifests  itself  as 
a  sort  of  relaxation  effect  due  to  the  incomplete  Legendre  expansion  attempting 
to  represent  the  sudden  ending  of  a  steep  temperature  rise.  The  amplitude  is 
small  enough  to  be  able  to  put  a  mean  curve  through  each  oscillation  and  the 
temperature  values  read  for  each  mean  would  only  vary  by  a  few  per  cent  from  those 
corresponding  to  the  oscillation.  Such  oscillations  are  typical  of  those  which 
can  occur  when  a  rapidly  changing  function  is  represented  by  the  first  few  terms 
of  its  Fourier  expansion.  It  is  possible  to  improve  the  convergence  of  Fourier 
series  and  hence  get  better  results  with  fewer  terms  by  using  the  so-called  sigma 
factors  to  modify  the  values  of  the  Fourier  coefficients.  If  necessary  these 
sigma  factors  could  be  used  to  smooth  out  oscillations  in  the  Legendre  representations 
and  this  would  provide  an  objective  means  for  introducing  the  mean  curves. 

Having  left  these  oscillations  without  any  smoothing  in  the  outer  surface 
temperature  profiles,  it  is  not  surprising  that  they  reappear  in  the  derived  heat 
transfer  profiles  as  shown  in  figures  24  to  29.  Here  again  the  amplitudes  are 
small  and  there  is  not  much  difficulty  in  guessing  at  positions  for  mean  curves. 

Outer  surface  heat  transfer,  Q  (b,t)  was  derived  from  the  time  histories  of  the 
Legendre  components  of  outer  surface  temperature,  l^fbjt),  through  use  of  equation 

(13).  Apart  from  the  heat  transfer  profiles  for  time  8  seconds  and  above,  where 
the  raw  data  is  noisy,  there  is  reasonable  agreement  between  the  profiles  over  the 
rear  part  of  the  hemisphere  for  the  three  meridian  planes  0=  *120°  and  0=0.  However, 
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over  the  front  part  of  the  hemisphere  there  are  marked  differences  in  the  heat 
transfer  profiles  shown  in  figures  24  to  29  and  the  question  is,  how  authentic 
are  these  differences?  This  question  is  important  in  deciding  where  boundary 
layer  transition  occurred.  The  initial  heat  transfer  profile  shapes  are  strongly 
influenced  by  the  levels  of  heat  transfer  at  0=0  and  0=20°  which  were  the  positions 
at  which  temperature  data  was  recorded.  The  curves  joining  the  points  0=0  and 
0=20°  are  Legendre  polynomial  interpolations  which  also  take  into  account  the  data 
recorded  from  all  the  other  thermocouples  in  the  same  row.  A  guide  to  the 
authenticity  of  the  derived  heat  transfer  differences  at  0=20°  between  the  three 
meridian  planes  0=  ±120°  and  0=0  can  be  obtained  by  comparing  heat  transfer 
variation  with  time  against  the  time-derivative  variations  for  the  smoothed  inner 
face  temperature  measurements.  The  temperature  profiles  of  figures  21,  22  and  23, 
show  only  slight  curvature  at  0=20°  and  hence  most  of  the  local  temperature  rise 
is  due  to  aerodynamic  heat  input  rather  than  circumferential  heat  conduction  and 
there  should  thus  be  a  correlation  between  the  time-derivative  of  the  inner  face 

temperature,  ^  (a,t)^  t^e  aerodynamic  heat  transfer  Qr  (b,t).  That  there  is 

such  a  correlation  can  be  seen  from  figure  30  and  this  means  that  the  differences 
in  the  derived  heat  transfer  results  for  0=20°  are  reflected  in  the  raw  temperature 
measurements.  The  scatter  in  the  derived  heat  transfer  results  for  the  stagnation 
point,  0=0,  is  shown  in  figure  31.  Three  different  results  for  stagnation  point 
heat  transfer  were  derived  because  each  row  of  thermocouples  was  treated 
independently.  The  scatter  can  thus  be  regarded  as  indicating  the  accuracy  of 
truncating  the  Legendre  polynomial  expansions  of  temperature  at  the  seventh  term. 

It  appears  from  figure  31,  that  the  markedly  different  profile  of  the  measured 
temperature  at  0=0,  0=20°,  has  carried  through  to  the  0=0,  0=0,  profile,  which  is 
seen  to  account  for  most  of  the  scatter  up  to  about  6.5  seconds.  In  the  light  of 
these  observations  we  can  now  re-examine  the  initial  heat  transfer  profiles  of 
figures  24  to  29  and  draw  some  conclusions  about  transition.  Up  to  time  4  seconds, 
figures  24,  25  and  26,  all  exhibit  a  characteristic  laminar  heating  profile  up  to 
about  0=15°.  At  time  5  seconds,  figures  24  and  26  show  profiles  for  0=0  and 
0=  -120°  which  indicate  transition  to  turbulent  flow  very  close  to  the  stagnation 
point.  This  situation  appears  to  exist  up  to  time  6  seconds,  if  the  possibility 
of  a  lower  value  of  Qr  (b,t)  than  that  shown  in  figure  26  for  0=  -120°,  0=0,  is 
accepted  in  accordance  with  the  amount  of  scatter  present  at  time  6  seconds,  shown 
in  figure  31.  For  0=  -120°,  figure  29  shows  a  strong  tendency  for  the  flow  to 
revert  to  transition  at  about  0=15°  for  time  6.4  seconds,  even  if  some  account  is 
taken  of  possible  scatter.  For  time  7  seconds  and  above,  it  appears  from  figures 
27,  28  and  29,  that  the  flow  becomes  fully  turbulent  at  about  0=20°. 

It  is  really  more  appropriate  to  look  at  Stanton  number  variations  when 
drawing  conclusions  about  transition.  Stanton  number  variations  are  shown  in 
figures  32  to  37.  These  figures  tend  to  support  the  conclusions  drawn  from  the 
heat  transfer  profiles  of  figures  24  to  29  and  the  initial  laminar  profiles  for 
times  less  than  8  seconds  in  figure  37  for  0=  -120°  are  more  pronounced  than  their 
heat  transfer  counterparts  in  figure  29.  The  two  Stanton  number  profiles  for 
times  of  8  seconds  are  greatly  distorted  versions  of  their  heat  transfer  counterparts 
because  at  these  times  the  difference  between  the  outer  surface  temperature  and 
recovery  temperature  rapidly  becomes  small  and  has  the  effect  of  magnifying 
oscillations  when  divided  into  the  heat  transfer  to  obtain  Stanton  number.  Also, 
the  Stanton  number  based  on  the  difference  between  recovery  temperature  and  outer 
surface  temperature  or  wall  temperature,  Tr  -  Tw,  stems  from  steady  flow  theory 
whereas  in  the  transient  case  with  a  falling  recovery  temperature  the  heat  transfer 
to  the  outer  surface  will  not  have  reached  zero  by  the  time  Tr  falls  to  Tw.  Hence, 
near  this  condition,  the  Stanton  number  tends  to  infinity.  Rough  calculations 
show  the  temperature  response  time  for  the  boundary  layer  well  back  on  the 
hemisphere  to  be  of  order  30  milliseconds  at  time  8  seconds.  A  comparison  between 
the  derived  Stanton  numbers  for  the  stagnation  point  and  those  for  the  steady 


WSRL-0G69-TR 


14 


laminar  flow  stagnation  point  theory  developed  from  reference  10  and  described  in 
reference  11  is  shown  in  figure  38.  In  figure  39  the  stagnation  point  heat  transfer 
calculated  from  this  laminar  theory  is  compared  with  the  derived  values  treating 
each  row  of  thermocouples  separately.  Agreement  is  seen  to  be  good  except  for 
about  one  and  a  half  seconds  past  maximum  heating  where  the  derived  values  are  up 
to  twenty  five  per  cent  higher  than  laminar  steady  flow  theory. 

Values  for  maximum  turbulent  heating  around  a  hemisphere  are  of  interest 
because  they  are  difficult  to  predict.  An  envelope  of  maximum  turbulent  Stanton 
number  values,  taken  from  the  heating  profiles  around  the  hemisphere  for  the  three 
meridian  planes,  is  shown  in  figure  40.  These  maxima  occur  at  about  30°  from 
stagnation  point  initially,  shifting  backwards  towards  60°  as  the  hemisphere 
becomes  hot.  In  this  experiment  the  maximum  turbulent  Stanton  numbers  ranged  from 
0.0022  to  0.0037. 


S.  CONCLUSIONS 

Outer  surface  temperature  and  heat  transfer  profiles  have  been  deduced  from 
inner  surface  measurements  of  the  temperature  in  a  hemispherical  shell  which 
reached  a  Mach  number  of  5.5  at  4  km  altitude  after  7  seconds  of  flight.  In 
deducing  outer  surface  conditions  a  new  mathematical  technique  has  been  tested 
and  found  to  have  been  effective.  In  any  future  experiment  using  this  technique 
it  is  recommended  that  at  least  nine  equally  spaced  thermocouples  per  row,  in 
addition  to  the  one  at  the  stagnation  point,  be  used,  and  that  finite  Legendre 
expansions  having  the  same  number  of  terms  as  the  number  of  thermocouples  in  a 
row  be  fitted  to  the  measurements  by  solving  the  resulting  simultaneous  algebraic 
equations  for  the  coefficients.  This  avoids  the  difficulty  of  having  to  perform 
numerical  integration  of  oscillatory  functions. 

The  derived  heat  transfer  profiles  are  characteristic  of  a  boundary  layer  flow 
which  is  laminar  in  the  stagnation  region  and  undergoes  transition  to  turbulent 
flow  within  the  first  20°  of  arc.  Agreement  with  laminar  stagnation  point  heat 
transfer  theory  for  steady  flow  is  good  except  for  just  after  maximum  heating 
where  the  results  from  the  experiment  are  up  to  some  25  per  cent  higher.  Maximum 
turbulent  heat  transfer  occurred  between  30^  and  60°  from  the  stagnation  point 
with  maximum  values  of  Stanton  number  ranging  from  0.0022  to  0.0037. 
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NOTATION 

inner  radius  of  hemispherical  shell 
outer  radius  of  hemispherical  shell 
thickness  of  a  flat  plate  conductor 
integers 

operational  transform  of  time  t 
radial  coordinate  in  shell  or  an  integer 
an  integer 
time 

distance  coordinate  in  a  flat  plate  conductor 

function  of  Spherical  Bessel  Functions  of  the  first  and  second  kinds 
defined  by  equation  (6) 

specific  heat  of  air  at  constant  pressure 

coefficient  of  thermal  conductivity  for  hemispherical  shell 
memory  coefficient  defined  by  equations  (14),  (15),  (16) 

Legendre  polynomial  of  order  2n 

rate  of  heat  flow  per  unit  area 

rate  of  heat  flow  per  unit  area  in  radial  direction 
Reynolds  number  based  on  shell  diameter 

Stanton  number  =  Qr  /  (Tf-  T^) 

temperature 
recovery  temperature 

wall  heated  surface  temperature 

free  stream  velocity 

argument  used  for  Spherical  Bessel  Functions 
thermal  di ffusivity 


S 


time  increment 


WSRL-0069-TR 


16  - 


6  angle  in  shell  coordinate  system  shown  in  figure  15 

M  cos  6 

f  1  dummy  time  variables 

p  air  density 

o  the  sth  zero  of  B 

s 

0  angle  of  meridian  plane  for  a  thermocouple 

array  as  shown  in  figure  6 

0  temperature  memory  coefficient  as  defined  by  equations  (8)  and  (9) 


Subscripts 


free  stream  conditions 
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Figure  1.  Photograph  of  rocket  test  vehicle  on  launcher 


(All  dimensions  in  mm) 


2.  Sketch  of  rocket  test  vehicle 
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Figure  3 


Figure  3.  Flight  path  showing  recovery  events 
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Figure  7.  Raw  temperature  measurements 


Figure  8.  Raw  temperature  measurements  {<t>=0,  0=20 
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Figure  15.  Coordinate  system  for  hemispherical  shell 
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Figure  16.  Amplitude  and  frequency  of  the  Fourier  conqponents  of  temperature 
(0=120°,  0=20°) 
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22.  Derived  outer  surface  temperature  distribution  (0=120  ) 
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Figure  24.  Derived  outer  surface  heat  transfer  distribution  (0=0,  t  <  6s). 
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Figure  26.  Derived  outer  surface  heat  transfer  distribution  (0=  -120°,  t  ^6s). 
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Figure  35.  Stanton  number  distribution  (0=0,  t  >6.4s). 
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Figure  38.  Comparison  of  derived  stagnation  point  Stanton  number  with 
laminar  theory. 
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Figure  39.  Comparison  of  derived  stagnat 
with  laminar  theory. 


Figure  40.  Peak  turbulent  heating  envelope. 
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